Exciting Dark Matter and the INTEGRAL/SPI 511 keV signal 
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We propose a WIMP candidate with an "excited state" 1-2 MeV above the ground state, which 
may be colhsionally excited and de-excites by e"'"e~ pair emission. By converting its kinetic energy 
into pairs, such a particle could produce a substantial fraction of the 511 keV line observed by 
INTEGRAL/SPI in the inner Milky Way. Only a small fraction of the WIMPs have sufficient 
energy to excite, and that fraction drops sharply with galactocentric radius, naturally yielding a 
radial cutoff, as observed. Even if the scattering probability in the inner kpc is <^ 1% per Hubble 
time, enough power is available to produce the ~ 3 x 10*'^ pairs per second observed in the Galactic 
bulge. We specify the parameters of a pseudo-Dirac fermion designed to explain the positron signal, 
and find that it annihilates chiefiy to e+e~ and freezes out with the correct relic density. We discuss 
possible observational consequences of this model. 
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I. INTRODUCTION 

One of the central unsolved problems of both par- 
ticle physics and cosmology is the nature of the dark 
matter. Although well motivated candidates exist, for 
instance a weakly-interacting massive particle (WIMP) 
or axion, they have thus far eluded definitive detection. 
This makes indirect signals, such as anomalous particle 
production in the Galaxy, especially interesting as a pos- 
sible indicator of the physics of dark matter. 



A. The Positron Annihilation Signal 

In 1970, the first detection of a gamma-ray line in the 
Galactic center found a line center at 473 ± 30keV, caus- 
ing the authors to doubt its true origin []|. Leventhal 
soon pointed out that the spectrum was consistent 
with positronium emission, because positronium contin- 
uum could confuse the line fit in a low-resolution spec- 
trograph, but worried that the implied annihilation of 
7 X 10^^ pairs/s was several orders of magnitude larger 
than his estimate of the pair creation rate from cosmic ray 
interactions with the interstellar medium (ISM). Since 
then, other balloon and space missions such as 
CGRO/OSSE \^] and most recently the SPI spectro- 
graph 0] aboard the International Gamma-Ray As- 
trophysics Laboratory (INTEGRAL) have greatly refined 
these measurements (see [1| for a summary of previ- 
ous measurements) . We now know that the line center is 
510.954 ± 0.075keV, consistent with the unshifted e+e" 
annihilation line but still find the estimated annihi- 
lation rate of ~ 3 x 10^^ e~^e~ pairs per second 0, Q 
surprising. The origin of these positrons remains a mys- 
tery. 

The OSSE data showed a bulge and a disk compo- 
nent, with ^ 25% of the flux in the disk[4]. The most 
recent INTEGRAL data show mainly the bulge compo- 
nent, approximated by a Gaussian with FWHM 6° 



The data are consistent with a disk component several 
times fainter than the bulge, when Galactic CO maps 
are used as a disk template. The spatial distribution of 
the ortho-Ps (positronium) 3-photon continuum and 511 
keV line are mutually consistent, with 92 ± 9% of the 
pairs annihilating through Ps The faintness of the 
disk component constrains the expected (conventional) 
production modes: cascades from cosmic ray interac- 
tion with the ISM, and nucleosynthetic processes during 
supernovae both would have a disk component. In or- 
der to achieve the observed bulge/disk ratio with type la 
SNe, one must assume that the positrons from the disk 
migrate to the bulge before they annihilate; otherwise 
the bulge signal would only be ~ 10% of that observed 

Indeed, Kalemci et al. 13" find that the quantity of 
pairs alone rules out SNe as the only source. Using SPI 
observations of SN 1006 (thought to be type la), they 
derive a 3f7 upper limit of 7.5% for the positron escape 
fraction. An escape fraction of 12% would be required to 
produce all of the positrons in the Galaxy, but even then 
they would be distributed throughout the disk, like the 
1.8 MeV line from ^^Al |lj|, rather than concentrated 



in the Galactic center. Hypernovae remain a possibility, 
if there are more than 0.02 SN 2003dh-hke events (SNe 
type Ic) per century [l^. GRB rates of one per 8 x lO'' 
yr X {Eg hb/W^^ erg) are also sufficient p^ . 

It is estimated that microquasars could contribute 
enough positrons to produce approximately 1/3 of the 



annihilation signal [17[, but again, it is not clear that 



their spatial distribution is so strongly biased towards 
the Galactic bulge. 

In summary, no obvious astrophysical explanation for 
this positron production exists; proposals in the litera- 
ture have difficulty explaining the size or spatial distri- 
bution of the signal, and often both. While it is still 
possible that conventional (or more exotic) astrophysical 
sources could produce the entire signal, the concentration 
of dark matter in the Galactic center motivates consid- 
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eration of an alternative hypothesis: that some property 
of dark matter is responsible for the signal. 

B. Dark Matter 

In recent years, WIMP annihilation scenarios have 
been invoked to explain a number of high energy astro- 
physical phenomena (see [T8| for a review) . It is possible 
to create the pairs directly from WIMP annihilation if 
the WIMP mass is a few MeV p^], but this mass range 
has no theoretical motivation. In extensions to the Stan- 
dard Model which addresses the hierarchy problem, such 
as supersymmetry, the WIMP mass is usually above the 
weak scale (100 GeV - 1 TeV) and the annihilation prod- 
ucts come out at energies very large compared to rUe- In 
the absence of a dense environment where a cascade can 
develop (column densities of ~ 10^^ cm~^ are required 
for the 7 — > e^e~ step), there is no way to partition the 
energy of these particles into the many thousands of pairs 
required for each WIMP annihilation. Consequently, it 
has generally been concluded that weak scale dark matter 
has nothing to do with the observed positronium annihi- 
lation signal. However, we shall see that WIMPs in this 
mass range could still play a role: they could convert 
WIMP kinetic energy into pairs via scattering. 

The kinetic energy of a 500 GeV WIMP moving at 
500 km/s is > 511 keV. If the WIMP has something 
analogous to an excited state - for whatever reason ~ in- 
elastic scattering could occur, raising one or both of the 
WIMPs to this excited state. If the energy splitting is 
more than 2me, the decay back to ground state would 
likely be accompanied by the emission of an e^e~ pair. 
One tantalizing feature of this mechanism is that it pro- 
vides access to the vast reservoir of kinetic energy in the 
WIMPs (roughly 10^" erg for the Milky Way), such that 
bleeding off 3 x 10''^ MeV/s = 5 x 10^^ erg/s reduces the 
kinetic energy by a negligible amount, even over the age 
of the Universe. Indeed, the unexplained excess in the 
bulge is much smaller than that. Therefore, this scat- 
tering would have negligible kinematic effect on the dark 
matter halo, except perhaps in its innermost parts. 

In this paper, we propose that the positrons in the 
center of the Galaxy arise precisely from this process of 
"exciting" dark matter (XDM). Remarkably, such a phe- 
nomenon is quite natural from the perspective of particle 
physics in theories with approximate symmetries in the 
dark sector, which are generally accompanied by nearly 
degenerate states, which can serve as excited states. Al- 
though the MSSM dark matter candidate, the neutralino, 
does not exhibit such a property, this is simply the result 
of the restricted field content of the MSSM, with no addi- 
tional approximate symmetries present. Other WIMPs, 
both arising in supersymmetry and from strong dynam- 
ics easily give this phenomenology ~ and usually without 
upsetting the desirable properties of WIMP models. 

The layout of this paper is as follows: in section |TT] 
we discuss the basic properties of such a scenario, and 



constraints on the model from the annihilation signal. 
In section Hill we review a simple example from particle 
physics with the required properties, utilizing a single 
pseudo-Dirac fermion as dark matter. In this scenario we 
will also comment on the annihilation rate into standard 
model particles and other direct and indirect tests of the 
model. 



II. CONSTRAINTS FROM THE 511 KEV 
SIGNAL 

While the spectral line and continuum data from SPI 
are robust when averaged over the entire inner Galaxy, 
the spectrum in each spatial pixel is noisy. In order to 
stabilize the fit and robustly determine the total positro- 
nium (Ps) annihilation rate in the bulge, the SPI maps 
are cross-correlated with spatial templates. Because the 
Ps emission in the inner Galaxy does not follow the stel- 
lar mass distribution, the choice of templates is somewhat 
arbitrary. Representing the disk with a CO map, and the 
bulge with a Gaussian, Weidenspointner et al. Q obtain 
continuum Ps fluxes in their "analysis bands" of fb = 
O-Se^ais X 10"^ cm-2 s"! and fco = ^■^'^t^^ x 
lO^'^ ph cm^'^ s^^ for the bulge and disk, respectively, 
finding 31% of the emission is in the bulge. Because the 
bulge component covers less solid angle than the disk, 
its surface brightness is much higher. Using the full SPI 
response matrix they then obtain a total Ps continuum 
(bulge plus disk) of (3.11 ± 0.56) x 10"^ ^m^^ 
and line flux of (9.35 ± 0.54) x lO"'^ ph cm^^ g-i. For 
a Galactocentric distance of 8 kpc = 2.4 x 10^^ cm, this 
implies a total pair annihilation rate of 1.1 x lO''^ s""'^, 
or 3.4 X 10'*^ s"^ for the bulge. Because of order 10% of 
that is explained by SNe [l^l, we take 3 x 10''^ s~^ as 
the excess pair creation rate in the Galactic bulge to be 
explained. 

SPI also measures the width of the annihilation line to 
be 2.37 ± 0.25 keV FWHM, which is consistent with a 
warm (7000 < < 4 x 10* K) weakly ionized medium 
p^ . Knowledge of the state of the ISM in the inner 
Galaxy may be used to constrain the injection energy 
of the e+ to less than a few MeV j20|, which agrees 
with the constraints set by Beacom et al. using internal 
bremsstrahlung [2l[ and propagation [13] ■ In our model, 
the few MeV constraint applies to the splitting between 
the states, not the mass of the WIMPs themselves. 

A. Form of the cross section 

In this section we use elementary considerations of the 
phase-space density and expected u-dependence of the 
cross section to derive constraints on the ratio of the 
splitting to the mass (S/M) from the slope of the 511 
keV emission. 

The velocity-averaged rate coefficient (i.e. number of 
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scatterings per time per density) (cru) is 

M(r)= I Cpvd^v'f[w,v)f{v',v)a{Vrel)Vrel (1) 



where /(v,r) is the phase-space density of particles at 
position r and velocity v, a{vrei) is the inelastic scat- 
tering cross section as a function of the velocity differ- 
ence Vrei ~ 1'^ ~ '^'L which vanishes below threshold 
{vrei < Vthresh), and the integrals are taken over all ve- 
locities. 

In order to make an estimate of the parameters in- 
volved, we approximate the velocity distribution as 
isotropic with no particle-particle velocity correlations 
(caused by e.g. net rotation of the DM halo). This re- 
duces {av) and / to functions of scalar r and v. The cross 
section of the proposed pseudo-Dirac fermion has a weak 
dependence on velocity, so we first focus on a simplified 
form of {av) where the velocity dependence is entirely 
encoded in the phase space of the final state particles, 
i.e., where the form of the cross section is 



aVrel 



(2) 



Here is independent of velocity and is the approx- 
imate cross section in the moderately relativistic limit 
(where the second factor goes to v). Here S is the split- 
ting of the excited state from the ground state and 
is the mass of the wimp. The threshold velocity for ex- 
citation is Vthresh = v^VMt^, SO the integral becomes 



{uv)(r) = ^^^cj,r.rj^v^f{v,r)dvj^v'^f{v\r)dv' 



2 

thresh 



with 



= + v'"^ + 2ot'cos( 



(3) 



(4) 



This conveniently separates the problem into a question 
of particle physics (cmr) and a question of astrophysics 
(the phase-space densities). 



B. DM velocity distribution 

The WIMP velocity distribution depends on the details 
of the formation history of the Galaxy, but for simplicity 
we use a Boltzmann approximation. Although there may 
be some WIMPs in the inner Galaxy moving faster than 
escape velocity, it is reasonable to assume that the dis- 
tribution goes to zero rapidly for v > Vesdf)- Therefore, 
we adopt 



N exp{-V^ /2vf^^) V<Vesc 
V > Vesc 



(5) 



with 1-D velocity dispersion w^ms = 200 km/s. N nor- 
malizes this so that J^"''" d^vf{v,r) = 1. Following Mer- 
ritt et al. [l^ we use a DM halo profile of the form 



Pq exp 



Da 



(6) 



where po 0.1 - 0.7 GeV [2J] is the DM density at the 
Solar circle (r — Rq), r_2 = 25 kpc is the radius at which 
the logarithmic slope of the profile is -2, and a « 0.2 is 
a parameter of the profile. This profile is inspired by a 
fitting formula for the logarithmic slope of density [2^: 



dlnp 
dlnr 



-2(r/r_2)" 



(7) 



and agrees well with simulations \2§\ . 

This profile only contains the dark matter, which does 
not dominate in the inner parts of the Galaxy. To esti- 
mate the escape velocity of the Milky Way, we note that 
the rotation curve is nearly flat with a circular velocity 
Vc ~ 220 km/s out to far beyond the solar circle [l^. In 
the Merritt profile the scale radius of r_2 contains much 
of the mass of the Galaxy, so we make the (conservative) 
approximation that the rotation curve is flat to 20 kpc 
and then there is no mass beyond that radius. A flat rota- 
tion curve implies that enclosed mass M(r) = rv'^/G, so 
p{r) ^ r~^. This density profile differs from the Merritt 
profile because it includes baryonic matter, which dom- 
inates the mass density in the inner part of the Milky 
Way. (The Black Hole dominates only inside the inner 1 
pc, so we neglect its influence.) 

The value of the one-dimensional velocity dispersion 
is chosen to agree with simulations |23| . Although some 
authors have correctly pointed out that in DM-only halos 
the dispersion decreases in the center (e.g. 28]), simu- 
lations including gas physics and star formation indicate 
an increase [13] . We defer a detailed discussion of these 
and related issues to a future paper. 



The gravitational potential $ 



-v^ at r = r_2, so 



$(r) 



therefore. 



GAf(r') 



dr' =t.2[ln(r/r_2)-l] (8) 



,(r) = V-2$(r) = 



Vc^2[l-\n{r/r^2)] 
Vcy/2r^2/r 



for r < r_2 
for r > r_2 
(9) 

431 km/s 



which for Vc = 220 km/s gives Wesc(8 kpc) 
and Vesc{^-5 kpc) = 674 km/s. Thus an energy thresh- 
old for the inelastic scattering corresponding to about 
i^rei = 1000 km/s will cause a sharp cutoff beyond a 
Galactocentric radius of 400 pc (about 3°) as observed 
by INTEGRAL. Other assumptions may be made about 
the mass profile outside of 20 kpc, but since any addi- 
tional velocities add in quadrature, they do not modify 
this result greatly. 
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FIG. 1: The fraction of pairs with velocity above Vthresh of 600, 800, or 1000 km/s, for Vrms ~ 200 km/s (solid) and Vrms ~ 
180 km/s (dashed). The approximate half-max radius of the observed 511 keV emission (vertical dotted) is shown for reference. 



It is instructive to consider the fraction of WIMP pairs 
with a relative velocity above threshold, as a function of 
radius. This provides a sense of how many pairs are 
participating in the inelastic scattering, independent of 
the detail of the particle physics model. This pair fraction 
is shown in Fig. [T] for three values of Vthresh and two 
values of Vrms, the 1-D velocity dispersion. 



C. Comparison with Observations 

To compare this model to the 511 keV map, we next 
integrate along the line of sight 



dl {a v){r) 71^ (r) 



(10) 



to obtain the number of pair creations per area per time. 
Figs, mo show the line of sight integral for various par- 
ticle model and halo profile parameters. For comparison, 
the most recent analysis of the INTEGRAL SPI signal 
finds a 6° FWHM Gaussian [ll| in agreement with pre- 
vious measurements by OSSE This corresponds to a 
pair intensity (inferred from both line and Ps continuum) 
of 



exp(-6'V2a2 



(11) 



where Lpair = 3 x 10^^ pairs/s, a = 0.044 rad (i.e., 
6°/2.355), i?o = 8 kpc, and 9 is the Galactocentric angle. 
This function peaks at / w 0.03 pairs s~^ cm~^ sr~^, and 
is shown in Figs. I2I3I for comparison. In these plots, we 
have attempted to employ standard values from a,po, 



and Vrms for the fiducial model. However, by taking 
more optimal values for all parameters, a cross section 
of (Tmr — 3 X 10~^* cm^ would be sufhcient to generate 
an acceptable rate to satisfy the INTEGRAL measure- 
ment. 

None of the model lines in Fig. [3] are in perfect agree- 
ment with the Gaussian, but keep in mind: 

• The Gaussian form chosen by the SPI team is only 
an approximation. A more thorough analysis would 
fit the SPI data directly to a class of XDM models 
and compare likelihoods, as in [28j . 

• The substantial instrumental backgrounds in SPI 
are removed by subtracting observations of "blank" 
sky 10° off the Galactic plane. We adjust the zero- 
point of the INTEGRAL measurements to match 
the model at 10°, improving the agreement (see 
upper dotted line. Figs. [2] and [3]). 

• We have chosen a narrow class of DM halo profiles 
for this benchmark model. A different profile could 
improve agreement with the Gaussian. 

• Significant variation of the velocity dependence of 
the cross section could arise in different particle 
physics models 

• Other factors could alter the results, e.g. the elastic 
scattering of WIMPs (see appendix) at some level 
may alter the entropy in the inner part of the halo 
and cause significant departures from the assumed 
velocity and radial distributions. 

Whatever deficiencies our model may suffer from in terms 
of the spatial distribution of the positrons, they are less 
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FIG. 2: Model curves for the constant sigma model in Eq taking amr = 1.3 x 10"^^ cm^ ^ 0.3 Gcv/cm j ^ 2.4 x 10"^'' - 

1.2 X 10~^^ cm^, both unsmoothed {thin) and smoothed by a 3° FWHM beam [thick) to approximate the spatial response of 
SPI. In all cases the solid red lines represent our fiducial model (M = 500 GeV, Vthresh = 850 km/s, and halo parameters Vrms = 
200 km/s, and Merritt index a = 0.2), and dashed blue lines represent variations of one parameter, (a) M — 400, 500, 600 GeV 
with S held fixed at 1 MeV so that Vthresh ~ 950, 850, 776 km/s, respectively, (b) vthresh = 800, 850, 900 km/s, keeping 
M = 500 GeV while 5 now varies above and below 1 MeV. (c) The Merritt profile index a is varied (see Eq[S|. (d) The RMS 
velocity in the inner Milky Way is varied. This is assumed to be constant with radius. The observed 6° FWHM signal (lower 
dotted line) and that signal plus an arbitrary baseline of 0.002 (upper dotted line) are shown in (a). The SPI zero is set by 
measurements 10° off the Galactic plane to correct for instrumental background. 



severe than annihilation models (e.g. |lj|]) because XDM 
provides a mechanism for a radial cutoff. Regardless of 
the details of f(r,v) this is an appealing aspect of XDM. 

The question remains whether such a model of dark 
matter could naturally arise from particle physics con- 
siderations. In the following section, we will see that 
such models occur simply in extensions of the standard 
model. 



III. MODELS OF EXCITING DARK MATTER 

It is quite straightforward to find models of this sort. A 
model of composite dark matter, for instance, involving 
a bound state of constituent particles, would be expected 
to have precisely this sort of structure, arising from its 
excited states |43j . However, one needn't invoke compli- 
cated dynamics in order to generate a splitting, as an 
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FIG. 3: Same as Fig. [2l but using the full expression for the inelastic scattering cross section (Eq. IA2|l derived in Appendix A, 
setting the coupling A+ = A_ = 0.18 and p = 0.4 GeV/cm^. 



approximate symmetry will suffice. 

An example of this is the neutron-proton mass differ- 
ence, which is protected by an approximate isospin sym- 
metry. Consequently, although the overall mass scale for 
the nucleons is GeV, the splitting is small, and protected 
against radiative corrections by the symmetry. We will 
take a similar approach here. 

Let us consider a Dirac fermion, composed of two Weyl 
fermions Xi,2- We will require the presence of a scalar 
field (p, which will be light, and mediate the process of 
excitation. The Lagrangian we consider is 

C = l^^<f>^^'cP + xU^^^'x^-mDXlX2 (12) 
-M(f>XlXl - - V{(p). 

This can be justified with a Z4 symmetry, under which 



Xi,2 e=^*''/2xi,2 and </> -</>. 

It is easiest to work in the basis x, X* — 1/V2(xi TX2) 
In this basis the mass matrix for the x's is 



M = 



(13) 



where A± = ^('^i i -^2)- There are simple conclusions 
to draw from this expression. At leading order, the 
Dirac fermion is understood as two degenerate Majorana 
fermions. If the Z4 symmetry is broken weakly to Z2, for 
instance by a expectation value, we expect these states 
to be split by a small amount 5 = 2A+ {(f)). It should 
be noted that because this is a symmetry breaking pa- 
rameter, it can be naturally much smaller than the mass 
of the dark matter particle. We shall consider a simple 



FIG. 5: Decay of the excited state into the ground state. 
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FIG. 4: Excitation diagrams for XX ^ X*X- 



Such a term induces a mixing Omix ~ ot {4>) /mh- Thus, 
decays of the excited state into the hghter state proceed 
through the diagram in Fig. [5] 



origin for the vev shortly, but for now will take it as a 
free parameter. 

The process of excitation is mediated by the Feynman 
diagram of Fig|4^. At very low velocity, one needs to sum 
the ladder diagrams |4j| . The calculation of the cross sec- 
tion is given in the appendix. Qualitatively, however, the 
cross section is often approximately geometric, i.e., set by 
the scale of the characteristic momentum transfer. That 
is, av ^ v/q^ ^ v/m5 ^ 10^^^ cm'^/s. This large excita- 
tion cross section, interestingly, does not occur for large 
coupling, but rather when the perturbative expansion in 
\^ /Att'^v breaks down. 

The origin of the splitting will come from an expecta- 
tion value of 0. This can be derived from a potential of 
the form 

m = + (14) 

from which we yield a vacuum expectation value (vev) 
^02^ _ m^/n. The mass splitting is then just S — 
2A+m0/^. 

The presence of such a vev seems to imply the existence 
of domain walls, which would be at the limit of what is 
allowed with regard to dominating the energy density of 
the universe [2^. This can be cured by the presence of 
cubic term a(j)^, or the promotion of the symmetry to a 
(continuous) gauge symmetry. Although the cubic term 
is aesthetically unappealing, as our focus here is on the 
phenomenology of such a model with regard to present 
astrophysics, we shall accept such a term and defer more 
appealing solutions to future work. 

At this point we have only explained the process of 
excitation, but have not included the production of e~^e~ 
pairs. This arises quite simply from a mixing term of the 
form 

a<l>^hh*. (15) 



A. XDM in the early universe 

While some worries, such as domain walls, can be cured 
simply, there are a number of non-trivial issues that must 
be addressed, most notably, the questions related to Big 
Bang Nucleosynthesis (BBN) and the relic abundance of 
the XDM. 

Let us begin by addressing the lifetimes of the various 
particles. Both excited and unexcited states of the dark 
matter will come into thermal equilibrium in the early 
universe, with the excited state decaying into the lighter 
state with an approximate lifetime 

V Hi (16) 

while the scalar (p decays with lifetime 

r^'^^ylsin'O^,^, (17) 

where is the Yukawa coupling of the heaviest fermion 
into which (j) can decay. (In general, as ^ IGeV, this 
is or y^^.) With reasonable values for the mixing pa- 
rameter, a scalar with mass m ^ IMeV will decay before 
nucleosynthesis. However, the decay of the excited state 
into the lighter state will naturally occur after nucleosyn- 
thesis, but well before recombination. Because of the low 
baryon to photon ratio, these positrons should thermal- 
ize before affecting the light element abundances, but 
it would be interesting to consider whether any effects 
would be observable. 

Although structurally somewhat different, the relic 
abundance calculation for XDM is ultimately very simi- 
lar to that of standard WIMPs. It is simplest, conceptu- 
ally, if one separates the sectors of the theory into a dark 
sector {x, X*; 0) ^^'^ the standard model sector. The rel- 
evant direct coupling between the two sectors is a. Nat- 
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FIG. 6: Diagrams contributing to thermal equilibrium in the 
dark sector and between the dark sector and the visible sector. 

uralness tells us that "^^/^^Lggs ^ 10^**, so as not to 
significantly correct the cj) vev. 

The relic abundance proceeds in a slight variation of 
usual WIMP freezeout. Here, x "^iH stay in equilibrium 
with (/), which, in turn, stays in thermal equilibrium with 
the standard model. Direct scatterings (fxj) hh are in 
equilibrium at temperatures above the Higgs mass for 
10~®, thus we can reasonably tolerate lOMeV 
from the naturalness constraint. In this way, the temper- 
ature of the dark sector is made equal to the temperature 
of the standard model particles. Whether this persists 
below the Higgs mass depends on a. The s-channel anni- 
hilation into standard model fermions shown in Fig. [5^ 
has a relativistic cross section aann = a'^m'j/{8TrmfJ. 
This should keep annihilations in thermal equilibrium 
down to T = TTip, nic ~ rrir, ruh for a ^ 10^^, 10^^, 10^^, 
respectively. 

The annihilation rate of x into occurs through the 
diagram in Fig [5Jd, with a cross section for xx ~^ 4'4> 
annihilation of (TannVrei = 8A^/(7rM^). The relic abun- 
dance from the freezeout of such an interaction is [29| 
(assuming equilibrium with the standard model) « 
2 X lO^^^GeV^^/ {(JannVrei)^ which yiclds an accept- 
able relic abundance for A ~ ^A//100TeV, or, roughly 
A ~ 1/10, which is precisely the region in which the ex- 
citation process is maximal. This result should not be 
surprising - it is the usual result that a weak scale par- 
ticle with perturbative coupling freezes out with the ap- 
propriate relic density. A rigorous study of the allowed 
parameter range is clearly warranted. 



B. Other decays and signals 

Although the freezeout process of XDM is similar to 
that of an ordinary WIMP, the scattering and annihila- 
tion products can clearly be different. One important 
feature is that the annihilation products are typically 
electrons and neutrinos for < 2m.„. However, unlike, 
e.g., MeV dark matter, when the particles annihilate the 
resulting electrons and positrons are extremely eiiergetic. 
This is intriguing because both the HEAT data W\ and 
the "haze" from the center of the galaxy [sil . Is^ point to 
new sources of multi-GeV electrons and positrons. Here, 
these high energy particles (from boosted on-shell <j) par- 
ticles) are related to the low energy positrons detected by 
INTEGRAL (from off-shell (j) particles). Such high en- 
ergy particles may create high-energy gamma rays from 
inverse scattering off starlight which could be observed 
in the future GLAST mission. We leave detailed analysis 
of these possibilities for future work. 

It should be noted that the electron and positron can 
be closed to allow a decay to two photons. Such a pro- 
cess is considerably suppressed, however, being down by 
a loop factor and suppressed by a. Consequently, these 
monochromatic photons should be produced at a rate 
roughly 10~^ that of the positron production rate, and 
are unlikely to be detectable within the diffuse back- 
ground (see, e.g., the discussion in [22j). 

IV. DISCUSSION - OTHER OBSERVABLE 
CONSEQUENCES 

A. Cluster heating 

In the most massive galaxy clusters, the velocity dis- 
persion is high (up to ^ 1000 km/s 1-D), and the major- 
ity of pair velocities may be above threshold. It is doubt- 
ful that the e^e~ signal could be seen by the current gen- 
eration of 7— ray telescopes, nor could the redshifted cos- 
mological signal be seen yet. However, the kinetic energy 
of the pairs (which is assumed to be < mg when they are 
created) must heat the intergalactic plasma. In the in- 
tergalactic medium (IGM), losses due to synchrotron and 
inverse- Compton scattering are negligible at these ener- 
gies. Losses are dominated by bremsstrahlung and per- 
haps Alfven wave excitation (see [1^), and both of these 
mechanisms heat the few x lO^K plasma in the IGM. 

The powerful X-ray emission from bremsstrahlung in 
galaxy clusters (e.g. Perseus, Virgo) must either be bal- 
anced by significant heating (perhaps from AGN) or else 
cooling matter must fall into the cores. Such cooling 
flows have been sought for years, but never found at the 
level expected. Recent X-ray spectroscopy has shown 
much less emission at 1/3 of the virial temperature 
than that expected from simple radiative cooling mod- 
els [s^l. Alternatives, such as light DM annihilation [sH] 
have been suggested, but once again XDM can provide 
similar results. Especially noteworthy would be models 



9 



of XDM with {av) rising as a higher power of w; these 
could produce a substantial fraction of the power needed 
to balance X-ray cooling in rich clusters. Furthermore, 
in such scenarios, the number of non-thermal (mildly rel- 
ativistic) electrons could distort the SZ effect in small 
but detectable ways. Because estimates of cluster heat- 
ing and SZ signal are strongly model-dependent, we defer 
discussion of this to a future paper. 



B. BH accretion 

The z ~ 6 quasars discovered by the SDSS [11] can 
be produced in simulations (Yuexing Li, priv. comm.) 
but only by seeding the simulations with 2OOM0 black 
holes at 2: ~ 30. A mechanism of dissipative scattering, 
such as XDM, provides the possibility of super-Eddington 
accretion. This could be important at high z, before 
baryons have had a chance to cool into potential wells, 
perhaps seeding the universe with black holes at early 
times. At later times, however, the more efficient cooling 



of baryons makes it unlikely that XDM accretion is a 
significant factor in BH growth in general. 



V. CONCLUSIONS 

We propose that dark matter is composed of WIMPs 
which can be coUisionally excited (XDM) and de-excite 
by e^e~ pair emission, and that this mechanism is re- 
sponsible for the majority of the positronium annihilation 
signal observed in the inner Milky Way. A simple model 
of a pseudo-Dirac fermion with this property is presented, 
though other models with similar phenomenology are 
likely possible. The XDM framework is strongly moti- 
vated by a specific problem, but has a rich set of implica- 
tions for cluster heating, BH formation, and other phe- 
nomena. The model naturally freezes out with approxi- 
mately the right thermal relic density, and annihilations 
today go mainly to e+e~. Although XDM is dissipative, 
it does not jeopardize stability of observed structure or 
disagree with any other current observations. 



APPENDIX A: EXCITATION CROSS SECTION 

The process of excitation arises from a sum of all ladder diagrams in addition to the tree level diagrams. We 
sum these using the spectator approach, where we take the intermediate fermion lines to be on-shell. Neglecting 
logarithmic dependencies on the momenta and scalar mass, one then finds the n*'' diagram to be proportional to 
q~^(A^/u)", where = _|_, A_A+ depending on the diagram considered. Here v is the velocity of one of the initial 
states X in the center of mass frame. This allows us to sum the diagrams by solving 



+ i- 



V 



(Al) 



and equivalently for the u-channel (crossed) diagrams. Here L is a loop factor which can depend on the states 
considered, and will be evaluated numerically. 

Because we are taking the intermediate fermion lines to be on-shell, we neglect diagrams where both fcrmions are 
excited, which are suppressed due to phase space. 



aVrel 



(A2) 



where arises from summing the internal single excitation ladder diagrams, excluding the external leg contributions 
to the amplitudes, which is 



(v^\^_\l ((Lo - Xi - 2X1 (Lo - L,f Xl +v' + XX (Lq - ^1)')) / 
'LlLfXl + 2XlLlLlXl + ((L?Ai + v^) LI + v^LI) A^ 
+2v'^X^_ {2Lq - Li) LiX\ +v^ + v^XiLj 



(A3) 



where here, Lo,i are the loop factors for diagrams involving zero or one excited internal fermions. f3 is the final state 
velocity of one of the dark matter particles. To a reasonable approximation, Li is just suppressed relative to Lq by 
phase space, that is, Li « LqP"^ /v^ . With this simplification 

M2 = v^X^_X\ + (t;2 _ /32)2 ^4 ^2 ^ (^2 _ /32)2 ^4 ^2 _ 2 (i,2 _ A^A^L^^ / (A4) 

{P^Ll {LlX\ + v^) Xt + 2p^X\Ll {2v^ - p^v^ + /J^Aj £2) + {L^XX + v^) {v^ + P*X\LI)) 
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In these expressions, we have kept terms to leading order in and used the non-relativistic approximation. 

We compute the factor Lq numerically, and find Lq ~ 3 ^ 10~^. Notice that in the perturbative limit, A'P — A^A^, 
as expected. This cross section peaks for Aj_ ~ Vrei/Lo when the lowest order analysis begins to break down, 
which is A ^ 10^^. For m = SOOGeV, A+ = A_ = 0.18 (the example used in the Fig. [3]) and v — 450 km/s 
(i.e., Vrei ~ 900 km/s), one finds avrei ^ .OlGcV^^. This is not a remarkable result, in that the exchanged is 
O(0.5GeV^), so in this region, the cross section is approximately geometric. I.e., one would estimate a cross section 
7r/0.5^ - the particles will scatter if the come within approximately d ~ 1/g. It is important to note, however, that this 
large cross section does not occur in the region where the theory is strongly coupled (i.e., A_ S> 1), but rather in the 
region where one needs to sum the perturbative series, /Attv ~ 1. Let us add that although the summation of the 
box diagrams does regulate the low-velocity properties of the scattering, when the cross section peaks (as a function 
of A), every term in the summation is equally important, and it is likely that a more detailed analysis, including all 
the relevant logarithms is necessary. 

Elastic Scattering Cross Section 

One can solve the equivalent equations to determine the elastic scattering cross section, as such a question may 
be relevant for questions of structure formation, as with the case of self-interacting dark matter. For the purposes of 
low-velocity elastic scattering, we can neglect excitation box diagrams, and one finds a total cross section 

a ^ l/ml. (A5) 

Given that can range from 1 MeV to 1 GeV we have cross sections from 10^^^ cm^ to 10^^^ cm^. With the very 
largest cross sections, the mean-free path of dark matter in the halo would be of the order of Mpc, which could be 
relevant for the questions of halo substructure [33|- However, these scatterings are at very low momentum transfer, 
with A9 ~ m^/M. Whether such small scatterings are significant enough to realize the scenario of Spergel and 
Steinhardt is worthy of additional study. 



Note Added: After an earlier draft of this work was completed, we became aware of concurrent work by M. Pospelov 
and A. Ritz [38)], in which a similar scenario was considered where the dark matter excites into charged states which 
decay via electron emission. In this model, excitation would occur via short-distance processes. Although the authors 
come to a negative conclusion as to the viability of such a scenario from those in this paper, we believe this is due to 
the lack of long distance excitations, and different assumptions about the halo model. 
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